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Abstract 

A Monte Carlo code, known as AASI, is developed for simulating energy spectra in 
alpha spectrometry. The code documented here is a comprehensive package where all 
the major processes affecting the spectrum are included. A unique feature of the code 
is its ability to take into account coincidences between the particles emitted from 
the source. Simulations and measurements highlight the importance of coincidences 
in high-resolution alpha spectrometry. To show the validity of the simulated results, 
comparisons with measurements and other simulation codes are presented. 
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1 Introduction 



Monte Carlo simulations have proven to be adequate tools for describing alpha- 
, beta- and gamma-particle transport, even in complex geometries. A great va- 
riety of computer codes have been developed for particle transport, dosimetry, 
particle physics and industrial applications. Different levels of sophistication 
exist among the codes, but even the simplest ones, which take into account 
Rutherford and Compton scattering, photoelectric absorption and continuous 
slowing-down of charged particles, can provide acceptable results. In many 
cases, simulation is the only practical way to explore the physics behind ob- 
served phenomena. 
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Alpha-particle spectrometry is a widely-used analytical method, for example 
in surveys of environmental radioactivity. The low activity of the samples ne- 
cessitates long counting times and a small sample-detector distance (SDD). 
The drawback of a small SDD is the possibility of coincidence summing be- 
tween the emitted alpha particle and subsequent emissions from the daughter 
nucleus. In addition, carefully designed sample preparation techniques are es- 
sential, since the alpha particles continuously lose their energy as they travel 
through matter. The energy loss leads to degradation of the spectrum quality 
via peak spreading, which increases with as the SDD is reduced. 

Simulations can be used to investigate the influence of various phenomena 
on the spectrum quality. The most important factors can be singled out and 
the measurement setup can be optimised. Moreover, unknown properties of 
the source, such as source density (or thickness) or source particle properties, 
can be determined. This is important, especially in the case of direct alpha 
spectrometry, when radiochemical sample treatment is omitted. The particle 
beam attenuation and interactions in basic research can also be examined. 

Many Monte Carlo simulation packages, such as the TRIM package [1], the 
GEANT software suite [2], and MCNP code [3], are suitable for simulating the 
alpha particle behaviour in the medium. However, these packages are not nec- 
essarily optimal for alpha spectrometry simulations. More specific approaches 
to alpha-spectroscopic simulations include the backscattering study of Ferrero 
et al. [4] and the investigation of aerosol particles by Pickering [5]. Roldan et 
al. [6] examined the spectrum quality at a small SDD. 

The present Monte Carlo simulation code, known as AASI (Advanced Alpha- 
spectrometric Simulation), is designed to simulate alpha-particle energy spec- 
tra. It is intended to be a comprehensive simulation package where all the 
major processes influencing the energy spectrum are included. Samples of var- 
ious types (aerosol particles, thick samples, non-uniform samples, etc.) are 
accommodated. Coincidences between the emitted particles are calculated us- 
ing nuclide-specific decay data that are stored in a library file prepared in 
extensible markup language, XML. Although the code has so far been applied 
to the simulation of alpha particle energy spectra from environmental samples, 
it can also be used for other applications. The typical running time on a 1.6 
GHz Pentium PC varies from seconds to a couple of minutes depending on 
the complexity of the simulation problem. The code is written in Fortran 95. 



2 Properties of the source and particle tracking 

Particle propagation through a material layer is determined by two physical 
processes: direction changes (scattering) and energy loss. The algorithm for 
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particle propagation in a given material layer proceeds as follows: 



(1) Emit a particle from a randomly selected point. 

(2) Calculate the distance, i.e. step length to the next scattering (or pho- 
toabsorption) event using the cross-section data. 

(3) If the particle is charged, adjust the step if a boundary of absorbing 
material is crossed. Calculate the continuous energy loss during the step. 

(4) If the particle energy is below the cut-off value, stop tracking. 

(5) If the particle crossed boundary of the material layer, proceed to the next 
layer if one is present. Otherwise, stop tracking. 

(6) Determine the next direction vector, i.e. scattering angles. 

(7) If the particle is a photon, determine the energy loss in the scattering or 
photoabsorption event. 



Here, characteristics of the source as well as the particle tracking method, i.e. 
determination of the scattering angles, are described. Calculation of the energy 
loss of alpha particles, electrons and photons is presented in the following 
sections. 



2. 1 Source 

Particle emission can originate from a point or from a finite-sized object. 
These objects, e.g. aerosol particles, can be embedded in the source matrix. 
The composition of the source matrix and the objects that emit radiation 
need not to be the same. For example, alpha particles can be emitted from an 
aerosol particle located inside a glass-fibre filter. In the spectrum simulations 
the number of alpha particle emissions is given in the input. 

The thickness of the source can be subjected to random fluctuation that is 
assumed to follow a Gaussian distribution with a user-given standard devia- 
tion a. To prevent impossibly large thicknesses, the resulting source thickness 
H a (r) is limited to 



where r is the radial position inside the source, s is a user-given parameter 
and H is the nominal (mean) thickness. 

Coordinates of the source particles are sampled as described by Siiskonen 
and Pollanen [7] . For sources with a random thickness and convex or concave 
sources, the vertical coordinate is sampled by the rejection method. Convex 
and concave source shapes are described with a paraboloid of revolution. The 
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user of the code supplies the central and side thicknesses of the source. If the 
source thickness is zero, all source particles lie on a plane. 

Source particles can have a spherical or elliptical shape. Spherical source par- 
ticles can have a log-normal size distribution. Inactive source particles can be 
coated with a uniform layer of radioactive material. In addition, a spherical 
shell of inactive material can be placed around a spherical source particle. 

The distance of the source particles from the source surface can be exponen- 
tially distributed inside the source matrix. This is a useful feature for investi- 
gating air filters in which radioactive aerosol particles are accumulated. This 
option is only available for cylindrical sources without thickness fluctuations. 
The distance d l of a source particle % from the source surface is obtained from 

cf = -Amf (2) 

where £ l is a random number between and 1, and A is the mean penetra- 
tion depth given by the user. Another user-given parameter, /, determines 
the fraction of the particles to be distributed according to Eq. (2). The rest, 
fraction 1 — /, is distributed on the source surface (d l = 0). Particles which 
have d l larger than the source thickness penetrate the source and are ignored. 
The total number of emissions from the source is reduced accordingly. 

An average solid angle subtended by the detector, the geometrical detection 
efficiency, is calculated. This is the number of hits received by the detector 
divided by the number of alpha particle emissions. The desired accuracy, the 
standard deviation of the efficiency, is given in the input. Calculation of the 
geometrical detection efficiency is necessary, for example, in direct alpha spec- 
trometry when radiochemical sample treatment is omitted. Tracers cannot 
then be used for quantitative activity determination. 

The measurement setup, consisting of the source, source backing, absorbing 
material layers and the detector, can be plotted in a file for visual inspection. 
Library routines for plotting were written by Kohler [8]. 

2.2 Particle tracking 

Electrons are tracked when they travel in the source, in the source backing and 
in the detector, including its dead layer. Photons are only tracked inside the 
detector (including its dead layer). Alpha particles are tracked in the source 
backing for backscattering studies. Otherwise, particles are assumed to travel 
in straight paths. The tracked particle is followed until it escapes the absorbing 
material or its energy falls below a cut-off value. When crossing a boundary 
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between two adjacent absorbing layers, the tracking step length is adjusted so 
that the step does not cross the layer boundary. 

Particle tracking starts with the sampling of the initial emission coordinates. 
The initial emission direction (0 O , <f>o) is chosen from a uniform distribution. 
Following the emission, the cosine of the polar angle 8 n (see Fig. 1) of the 
tracked particle is determined by 

cos 6 n = cos uo cos 6 n -i + sin uj sin 6 n -i cos ip (3) 



where 9 n is the polar angle after nth scattering, u is the scattering polar angle 
and ip is the scattering azimuthal angle. The cosine and sine of the azimuthal 
angle <p n are given by 



sin0 n = sin(0 n - 4> n -i) cos0„_i + cos(0 n - (f) n -i) sin0 n _i 
cos n = cos(0 n - n _i) cos 0„_i - sin(0„ - 0„_i) sin n _ x 



(4) 
(5) 



where 



■ (a a ^ sinwsin^ 

sm(0 n - 0„_i) = — (6) 



'n 



and 



COS(0„ 



cos u — cos 6 n cos 6. 



n-1 



'n-1 — 



sin 6 n sin 6 n - 



(7) 



All angles are in laboratory coordinates. The scattering angle u depends on 
the differential scattering cross section. After the initial emission, the particles 
undergo successive scatterings which are assumed to be statistically indepen- 
dent. 



Alpha particle scattering is calculated in the centre-of-mass frame. Before 
the determination of cos# n , the scattering angle is transformed to laboratory 
coordinates via 

siny? 

tancu = — — (8) 



where (p is the scattering polar angle in centre-of-mass coordinates, m a and 
M are the masses of alpha particle and target atom, respectively. 
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3 Simulation of alpha particle behaviour in medium 



3.1 Energy loss 

Alpha particle energy loss is calculated as described in Ref. [7], using the 
stopping power parametrisation of Ziegler as described in Ref. [9]. The total 
stopping power is the sum of the stopping power due to electrons, S e , and the 
nuclear stopping power, S n . In the energy region of interest (below 10 MeV), 
S e is parametrised as 

111 

»Jc »Jlow »Jhigh 

where 

5 low = Cl E^ + c 3 ^ 4 (10) 

and 

Shigh = c 5 Ep In (c 7 E~ 1 + c 8 E p ) . (11) 

Values of parameters ci, . . . , c$ are tabulated in Ref. [9]. Here, E p is the energy 
of a proton moving at the same velocity as the alpha particle in question. For 
composite materials, other parametrisations are also available [7]. 

An arbitrary number of absorbing material layers can be added between the 
source and the detector. The user supplies the number of layers, their atomic 
and mass numbers, densities, thicknesses and standard deviations representing 
the thickness fluctuations. Alpha particles are assumed to travel straight paths, 
except in the source backing where their path is tracked collision by collision. 

3.2 Energy loss straggling 

Straggling of the alpha particle energy loss can be approximated by a Gaussian 
energy distribution. Although not strictly correct with thin absorbing layers 2 , 
it gives a reasonable estimate in many cases. Standard deviation o"g of the 
Gaussian distribution depends on the maximum energy transfer in one collision 
between the alpha particle and an atomic electron, E max , given approximately 

2 The asymmteric Vavilov (or Landau) distribution is better with thin absorbers. 
However, computationally they are much more complex. 
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by 2m c c 2 f3 2r ) 2 . The parameter 7 is defined through r jm a c 2 = E, where E is 
the energy of the alpha particle. The deviation is given by [2] 

ol = E max S(l - y) (12) 



where (3 is the alpha particle velocity in units of c. Parameter 5 is the average 
energy loss in the material layer in question, 

5 = 0.0614-f-p 5x keV. (13) 
p 2 A 



Here, Z and A are the atomic and mass number of the target, respectively, p is 
the material density in g/ cm 3 and Sx is the distance travelled in micrometers. 



3.3 Scattering in the source backing plate 



Alpha particles can be tracked in the source backing plate. Screened elastic 
Rutherford scattering is used to determine the changes in the flight direction. 
Between the elastic collisions, alpha particles are assumed to lose their energy 
continuously. The mean free distance between the collisions is calculated from 
the potential 

V(r) = JL-t±**l\ (14) 
47re r 

where a is the screening radius, e is electron charge magnitude, e is the 
permittivity of free space and r is the radial distance. The resulting total 
cross section is 

Z 2 

o~(E) = 7i (hca) 2 — - — (15) 



where a is the fine structure constant and screening parameter 

V ={l)s^E- (16) 



The screening radius is given by [10] 
0.885a 



4 3 + ^ 2/3 



(17) 
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where ao is the Bohr radius. The effective charge of the alpha particle, z e g, is 
calculated as described in [9,11]. The mean free distance (step length) between 
the collisions, L, is sampled from 

L = -(a(E)Ny l \nt, (18) 



where N = N A p/A is the atomic density and N A is the Avogadro constant. The 
mean atomic spacing iV~ 1//3 is used as a step length if L < iV -1 / 3 . Moreover, 
if the energy loss between two successive collisions is more than five percent 
of the alpha particle energy, the step length for the energy loss calculation is 
reduced until the loss is less than five percent. 

Angular deflection in the scattering event from a potential (14) is given by 
coscj = 1 ^ — (19) 

and ip = 27r£ (the £'s are independent random numbers). 



3.4 Detector response to alpha particles 



Alpha particle energy loss in the detector dead layer is treated as described 
in section 3.1. When the alpha particle hits the active volume of the detector, 
all its remaining energy is assumed to be deposited. In other words, alpha 
particles are neither tracked nor is their energy loss calculated in the active 
volume of the detector. Instead, a simplified solution is chosen which notably 
reduces the calculation time. 

The properties of the detector are read from the user-prepared file. The param- 
eters are the atomic and mass numbers of the detector material, detector ra- 
dius and thickness, dead layer thickness, detector full-width at half-maximum 
(FWHM) and the parameters of the exponential tailing function. 

Measurements show that the detector response to monoenergetic alpha parti- 
cles is not Gaussian [12,13,14]. To take this into account, a double-exponential 
tailing function can be added to the detector response. The resulting energy 
E is sampled from a distribution 

P( E ) = -J^L e -v(Eo-E) + _V_ e -„(E -E)^ E < Eo (2Q) 
1 + R 1 + R 



where E is the energy of the incoming alpha particle. The user supplies the 
parameters z/, \i and R. They should be determined from the measurements of 
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good quality (i.e. thin) sources at a large SDD. Parameter R > is the ratio 
of the areas of the two exponential distributions. Typical values for Canberra 
PIPS with an area of 450 mm 2 are v — 0.1 keV -1 , \x = 0.02 keV -1 and 
R = 12.0. FWHM of the Gaussian detector response is 14 keV. Convolution 
with the Gaussian detector response is done after the tailing. 



4 Simulation of electron and photon behaviour in medium 



4-1 Electrons 



Large-angle deflections of the electrons result from the screened elastic Ruther- 
ford scattering, Eqs. (14) and (15), with Z replaced with Z(Z + 1)/2, see, e.g. 
Ref. [15]. The screening parameter 77 can be chosen from three alternative 
models. Nigam et al. [16] suggested that 



?7n = 5.448- 



E P 



(21) 



where E e is the electron kinetic energy. Adesida et al. [17] fitted electron 
scattering data in aluminium and proposed that 



VA = 2.61- 



£2/3 



(22) 



Moliere [18] (see also [16]) concluded that 



Vm = a 



v2/3 0.36 + 1.20(aZ//3) 2 
r(r + 2) 



(23) 



where r = E c / (m c c 2 ) and m e is the mass of the electron. 

Between the elastic collisions the electrons continuously lose their energy. The 
energy loss is calculated by using Bethe's formula [19] 



dE e 
dx 



' e 2 \ 2nN A Zp 
4^J m e c 2 (3 2 A 



E c (E c + m e c 2 ) 2 p 2 
2I 2 m P c 2 



In 



+ 1 - P 2 - In 2 (2y/l -P 2 -l+[3 2 )+ \(1 - y/l - (3 2 ) 



(24) 



where / = 16Z 9 (eV) is the average ionisation energy of the target atom [2]. 
We neglect bremsstrahlung and the production of delta electrons and X-rays, 
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since their influence on the detected electron energy spectrum is small. The 
mean distance between the collisions and the angular deflection is calculated 
as for the alpha particles, except that the maximum energy loss per step is 
ten percent of the electron energy. 

When an electron hits the detector, its path is followed through the dead layer 
and into the active volume. In the dead layer, the user has an option to have 
a partially-depleted region, where the electron deposits part of its energy in 
the detected signal. The amount of energy deposited in the formation of the 
signal, E s , is then given by 

E s = AE C -^—, 0<z< 5z dl (25) 

0Z d \ 

where AE e is the electron energy loss at depth z and 8z A \ is the dead layer 
thickness. In the active volume of the detector, the deposited energy is equal 
to the detected signal, E s = AE e . 

The number of backscattered and transmitted particles from the detector is 
calculated. A particle is counted as backscattered when it escapes from the 
front side of the detector, and as transmitted when it escapes from other sides 
of the detector. For backscattering studies, a parallel electron beam hitting 
the detector surface perpendicularly can be used. 

4.2 Photons 

Photons are assumed to interact via photoelectric absorption and Compton 
scattering. Pair production is ignored, since we are interested in low-energy 
phenomena. The mean free distance is calculated from the total cross section 
of the above-mentioned interactions. A random number is used to decide which 
interaction occurs at the interaction point. 

The total photoelectric absorption cross section is read from a text file and 
interpolated. Data were obtained from the National Institute of Standards and 
Technology database [20]. If data for the element in question do not exist, an 
analytical approximation for the cross section [21] 

MEJ = (j^j 1 ' 'al (26) 

is used. Here, E\ = 0.0136 keV and E 1 is the photon energy in keV. This 
approximation overestimates the total cross section for most elements, espe- 
cially at low energies (less than 100 keV). For example, the overestimation is 
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approximately a factor of 10 in Si when E 1 is between 15 keV and 100 keV. As 
the overestimation is quite large, the user should supplement the photoelectric 
data library for the element in question, if possible. After the photoelectric 
absorption, an electron is ejected in the direction of the incoming photon, with 
E e = E y . X-rays produced in this process are ignored. Their energy is often 
so small that the photon is absorbed in the detector. 

The differential cross section for Compton scattering is 



dac (ahc) 2 1 

~dn~ 2 [m c c 2 + £ 7 (1 -cos to)} 2 

x ^\ -— - - + l + cos 2 cu (27) 

[ m c c z [m c c z + Ej{1 — cos lo)\ J 

neglecting the binding energy and momentum of the atomic electron. The 
photon scattering angle is calculated from distribution (27) using the rejection 
method by Brusa et al. [22]. The total cross section for Compton scattering is 



m e E^ 



2 / 2 \ 2 

m c c / m c c \ 
1 -2— 2 ' 



E, 



(m c c 



7 . 
2\2 



In 



2E^ + m c c 2 
m c c 2 



1 ^rn e c 
+ 2 + "E^ ~~ 2(2£ 7 + m e c 2 ) 2 



(28) 



The total interaction cross section is <j(E y ) = a P (E y ) + a c (E 1 ). The mean 
free distance is then 



L 1 {E 1 ) = -{a{E^N)- 1 \ni. 



(29) 



5 Treatment of coincidences 



The majority of alpha emitters have a significant decay branch to excited states 
of the daughter nucleus. The excited states decay by gamma-ray or conversion 
electron emission. Since the lifetimes of the excited states are typically much 
shorter than the integration time of the data acquisition electronics, pulse 
summation between the alpha particle and particles emitted by the daughter 
nucleus may occur. The summation is more pronounced when the SDD is 
small. 

The summation may lead to distortion of the peak shape and, thus, may have 
an influence on nuclide identification. A good example is separation of 239 Pu 
and 240 Pu, which is difficult even in the case of a high-resolution detector 



11 



and a sophisticated spectrum deconvolution code. Another example is the 
coincidence summation of 241 Am, whose main alpha decay branch leads to 
third excited state of the daughter 237 Np, resulting in a clearly visible bump 
on the high-energy side of the 241 Am main alpha peak. 

The probability of each alpha decay branch is given in the nuclide library file, 
consisting of a schema file and actual library in XML format. The Fortran- 
XML interface is written by Markus [23]. The alpha decay branch for an 
individual decay is selected using a random number. The nuclide library file 
also contains decay routes of the excited states of the daughter nucleus. Each 
decay route has a known probability of occurence (yield), decay type (gamma 
or conversion electron emission), initial and final state indices, and energy of 
the emitted particle. The emission of a conversion electron is associated with 
an X-ray, whose energy is given in the library. 

For each excited state of the decay route a random number is used to select 
the next decay channel (i.e., final state and emitted particle). The route is 
followed until the ground state is reached. For each emitted particle, the emis- 
sion direction is sampled. After the conversion electron emission, an X-ray is 
emitted before the cascade is followed further. We assume that each conversion 
electron is associated with X-ray emission. This is a simplification, since we 
overlook fluorescence yields and Auger electrons. The approximation is good 
for heavy elements, whose K-shell fluorescence yields are close to 100%. 

If the particle deposits energy in the active volume of the detector, a coinci- 
dence is formed and deposited energy is added to the alpha particle energy. If 
cascade consists of n subsequent decays, the alpha particle can be in coinci- 
dence with m < n particles. Deposited energies of those m particles are then 
added to the alpha particle energy. 

The algorithm to calculate the coincidences proceeds as follows: 

(1) Check that decay routes exist, i.e., transitions are available for the present 
state. If no route is found, exit the loop. 

(2) Use a random number to select the decay route, i.e., decay type, line 
energy and final state index. 

(3) If the emitted particle is an electron, follow it through the source and its 
backing. Determine if the particle travels towards the detector. 

(4) If the particle hits the detector, simulate the deposited energy. 

(5) If particle deposits energy to the detector, a coincidence is formed. Add 
the deposited energy to the alpha particle energy. If the particle is backscat- 
tered or transmitted, increase corresponding counters. 

(6) If the emitted particle is an electron, emit the associated X-ray. Go to 
(3). 

(7) Go to (1). 
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The lifetimes of the excited states are available in the library file. However, 
when coincidences are calculated, they are not taken into account. The life- 
times are assumed to be short enough, compared to the integration time of 
data acquisition electronics, for a coincidence to be seen. 



6 Comparisons with measurements and other simulations 

Geometrical detection efficiency and the alpha particle energy loss were inves- 
tigated by Siiskonen and Pollanen [7]. They found an excellent agreement with 
earlier results and measurements. To further confirm the homogeneity of the 
emission point distribution inside a source, we compared a simulated alpha 
particle energy spectrum from a thick sample with one obtained by numerical 
integration. In the comparison, a parallel alpha particle beam was considered 
(corresponding to a very large SDD, polar angle 6 = 0) in order to keep the 
numerical integration tractable. The agreement between the simulated spec- 
trum and the one from numerical integration is good (Fig. 2). This confirms 
the homogeneity of the emission point sampling, which is qualitatively also 
shown in Ref. [24]. Equally good agreement was obtained when the source 
was assumed to be spherical in shape (results are not shown here). 

Electron backscattering can be used to examine the quality of electron trans- 
port, since backscattering is sensitively depent on continuous energy loss and 
angular deflections in elastic Rutherford scattering events. Electron backscat- 
tering coefficients for various elements are compared in Table 1. The agreement 
between the experimental values and those of the present work is good. How- 
ever, the results of the present work are higher than the simulated results of 
Gueorguiev et al. [25]. The difference could be explained by different calcu- 
lation of the continuous energy loss. Gueorguiev et al. used different average 
ionisation energy and a three-point difference scheme. The convergence of the 
present backscattering calculations was ensured to two significant digits. 

To investigate the influence of coincidences, measured spectrum of 241 Am is 
compared to the simulation. The coincidence summing of alpha particles with 
photons and electrons from 237 Np is clearly visible as a bump above 5490 keV 
(Fig. 3). This bump is absent at a large SDD. Since experimental subshell 
conversion coefficients for 237 Np were not available, the relative yields for con- 
version electrons were set as follows: the yield of the L a line was assumed 
to be 20% and that of was assumed to be 80%. Figure 3 illustrates that 
simulations are accurately able to explain the effect of coincidence summing. 
When the coincidences are ignored in the simulation, the resulting spectrum 
clearly disagrees with the measurement (Fig. 4). 
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7 Discussion 



The present Monte Carlo simulation code, known as AASI, is designed for 
simulating energy spectra in alpha spectrometry. The code was originally de- 
veloped for estimating the influence of source characteristics on the alpha par- 
ticle spectra, for example in the case when the source quality is not optimal for 
high-resolution alpha spectrometry. The sources may be considerably thicker 
than those obtained from radiochemical sample treatment and their thickness 
may not be uniform. The source may even be an aerosol filter in which radioac- 
tive materials are deposited. This option is useful when the presence of alpha 
particle emitting materials in the filter must be identified rapidly, i.e. alpha 
particles are counted directly without prior radiochemical sample manipu- 
lation [7]. This information may be of utmost importance should a nuclear 
incident or malicious dispersal of radioactive material into the environment 
occur. 

Later development of the code is focused on the development of modelling of 
the detector and the effect of alpha-electron and alpha-photon coincidences 
on the measured spectra. The code can be easily used for various alpha de- 
tectors provided that the detector response can be treated as the convolution 
of a Gaussian part of a peak and a double-exponential low-energy tail (Eq. 
20). The comparison of measured and simulated spectra, especially in the case 
of a thin 241 Am source, highlights the importance of coincidence phenomena. 
When the source-detector distance is small (less than approximately a few cm) 
the coincidences must be taken into account in unfolding the spectra. Quanti- 
tative separation of nuclides such as 239 Pu and 240 Pu in alpha spectrometry is 
questionable if coincidences are neglected. Developments in the simulation of 
electrons and photons in the medium also facilitate the use of AASI in basic 
research. The detector response to various particle beams can be examined. 

Future development of the code will be addressed in the construction of an 
appropriate library data file and a user-friendly interface. The validity of the 
results, i.e. the agreement of simulated and those obtained from other calcula- 
tions or measurements, is verified here and in previous publications. However, 
validation is a continuous process not limited to this article only. 

An executable version of the code, the input files and a library file needed for 
the execution of the code are available from the authors upon request. 
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Table 1 

Comparison between experimental (Exp.) and simulated (MC) electron backscat- 
tering coefficients for various elements when E e = 30 keV and for normal incidence. 
In all simulations, the screening model of Nigam et al. [16] was used, see Eq. (21). 



Element 


C 


Al 


Cr 


Au 


Exp. [26] 




0.144 




0.482 


Exp. [27] 


0.060 


0.155 


0.270 


0.521 


MC [25] 


0.042 


0.124 


0.231 


0.502 


MC, this work 


0.053 


0.142 


0.253 


0.511 
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A 



Fig. 1. Coordinate system. The original particle direction is OA, the new direction 
is OB. The scattering polar angle is u and the azimuthal angle is tp. Laboratory 
polar angles are 9 n and 9 n -±. The z-axis is parallel to the detector symmetry axis 
and points from the source towards the detector. 
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Fig. 2. Comparison between simulated alpha particle energy spectrum from a convex 
source (side thickness ^m, central thickness 2 fim, solid black line) with a spectrum 
obtained by numerical integration (dashed grey line) from the same source. Alpha 
particles were assumed to travel in parallel tracks. 
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Fig. 3. Measured (dots) and simulated (solid line) spectrum of 241 Am. The 
source-detector distance was 5 mm and a detector FWHM of 14 keV was used 
in the simulation. Detector tailing parameters were v = 0.1 keV -1 , n = 0.02 keV" 1 
and R = 12.0. 
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Fig. 4. Simulated spectra of 241 Am when coincidences are taken into account (solid 
line, same as in Fig. 3) and when coincidences are ignored (circles). See the caption 
of Fig. 3 for simulation parameters. 
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